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ABSTRACT 

Context. Small inner working angle coronagraphs, like the vortex phase mask, are essential to exploit the full potential of ground- 
based telescopes in the context of exoplanet detection and characterization. However, the drawback of this attractive feature is a high 
sensitivity to pointing errors, which degrades the performance of the coronagraph. 

Aims. We propose a tip-tilt retrieval technique based on the analysis of the final coronagraphic image, hereafter called Quadrant 
Analysis of Coronagraphic Images for Tip-tilt Sensing (QACITS). 

Methods. Under the assumption of small phase aberrations, we show that the behaviour of the vortex phase mask can be simply 
described from the entrance pupil to the Lyot stop plane by Zernike polynomials. This convenient formalism is used to establish the 
theoretical basis of the QACITS technique. Simulations have been performed to demonstrate the validity and limits of the technique, 
including the case of a centrally obstructed pupil. 

Results. The QACITS technique principle is validated by experimental results in the case of an unobstructed circular aperture, and 
by simulations in presence of a central obstruction. The typical configuration of the Keck telescope (24% central obstruction) has 
been simulated with additional high order aberrations. In these conditions, our simulations show that the QACITS technique is still 
adapted to centrally obstructed pupils and performs tip-tilt retrieval with a precision of 5 x 10~^A/D when wavefront errors amount to 
/l/14rms and 10~^A/D for 2/70 rms errors (with A the wavelength and D the pupil diameter). 

Conclusions. We have developed and demonstrated a tip-tilt sensing technique for vortex coronagraphs. The implementation of the 
QACITS technique is based on the analysis of the scientific image and does not require any modification of the original setup. Current 
facilities equipped with a vortex phase mask can thus directly benefit from this technique to improve the contrast performance close 
to the axis. 

Key words. Techniques: high angular resolution. Methods: analytical, Methods: numerical 


1. Introduction 


Vortex coronagraph s (VC, iMawet et al.ll2005l: lEoo et alll2005l: 
Mawet et M l2009l) stand amongst the most promising focal 
plane phase masks envisioned for the next gene ration instru¬ 
ments of future very large telescopes (e.g. METIS, iBrandl et~^ 
2Q14I) . Theoretically, this coronagraphic solution provides a 
perfect star light rejection and other valuable features for di¬ 
rect imaging and characterization of exoplanets: achromaticity, 
continuous 360° discovery space and small inner working an¬ 
gle (angular distance where the off-axis transmission reaches 
50%). Eor these reasons, vortex phase masks already equip 
several infrared instruments on 10 m class groun d-based tele- 
scopes, namely VLT/NACO (iMawet et"^ l2013h . VLT/VISIR 
( Delacroix et all l2()12l: iKerber et all l2()14l). LBT/LMIRCam 
( Defrere et all l2014l) . Subaru/SCExAQ~ iljovanovic et al.l l201^ 


and very recently Keck/NIRC2. Scientific results have been ob¬ 
tained using the coronagraphic mode of these facilities, lead- 
ing to the detection of exoplanets and circumstellar disk s (e.g. 
lAbsil et aP 120131: iMilli et all l20l4 iReggiani et all l2014l) . The 
off-axis well-corrected subaperture on the Palomar Hale tele- 
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scope (ISerabvn et al.ll2007l) also provides a VC mode, which has 
led to impressive results (ISerabvn et n]l2010l) . including the de¬ 
tection of a companion very close to its host sty (e C ephei), at 
an angular separation of 1.12/D (iMawetet al .1120111 2 and D 
being the wavelength of observation and the telescope diameter 
respectively). 

However, a small inner working angle comes inevitably at 
a cost. Vortex phase masks, and in particular vortices of topo¬ 
logical c harge In = 2 like the Annular Groove Phase Masks 
(AGPM, IMawet et~^l2005l) . are amongst the focal plane masks 
that offer the narrowest inner working angle (down to 1 2/D), 
but it also means that they are highly sensitive to the centering 
of the star on the mask. Accurate tracking systems are there¬ 
fore required to limit the contrast loss due to pointing errors. A 
variety of low-order aberration sensing tec hniques exist s and i s 
used in current instruments, as reviewed bv IMawet et al.l (l2012h . 
In order to avoid non-common path errors, the sensor must be 
placed as close as possible to the coronagraphic phase mask. 
Solutions include sensors built just before the coronagraphic 
mask, like the Diffe rential Tip-Tilt Sensor (DTTS) of SPHERE 
(iBaudoz et al.ll2Qld where part of the light is diverted thanks to 
a dichroic plate), or the Cal low-order wavefront sensor of GPI 
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Entrance pupil Focal plane Lyot plane 



Fig. 1. Standard coronagraph layout, with a vortex phase mask at the 
focal plane and a Lyot stop at the second pupil plane (Lyot plane). 


([Wallace et al.ll201Q[ which uses the light rejected by the cen¬ 
tral spot of the occulting mask). For focal-plane phase masks, 
the later solution cannot be i mplemented , but a comparable so¬ 
lution has been proposed bv ISingh et al.l (l2014l) . making use of 
the light rejected by the coronagraph thanks to a reflective Lyot 
stop. Finally, phase retrieval techniques can also be applied di¬ 
rectly from the image acquired by the scientiflc detector, like 
the COFFEE sensor implemented in the SPHERE instrument 
(ISauvage et al.ll20l3) . 

In this paper, we propose a solution belonging to the latter 
category. It is based on the analysis of the flnal image produced 
by a VC to retrieve the tip-tilt affecting the beam incident on 
the phase mask. The principle of this technique, referred to as 
Quadrant Analysis of Coronagraphic Ima ges for Tip-tilt S ens- 
ing (QACITS) has first been introduced b y [Mas et al.l jl2Q12l) for 
the four quadrant phase mask (FQPM, iRouan et M l200Qh . It 
consists of quantifying the asymmetry observed in the corona- 
graphic point spread function (PSF), using the same principle as 
a quadrant cell detector. The differential intensities, or intensity 
gradients, are related to the pointing error and allow the estima¬ 
tion of the tip-tilt aberration affecting the beam. The simplicity 
of this technique makes it very easy to implement on current in¬ 
struments working with a vortex phase mask, as there is no need 
for any modification of the optical setup. 

In the next section, we describe the QACITS technique ap¬ 
plied to the perfect VC with an unobstructed pupil, and in partic¬ 
ular the mathematical model linking the asymmetry in the image 
and the tip-tilt, followed in Sect.[3]by an experimental validation 
of the model. For the sake of clarity, the details of the analytical 
computation are given in the appendices, where we introduce a 
formalism based on Zernike polynomials. In Sect.lH we detail 
the implications of a central obstruction on the PSF shape, and 
thus on the model used in QACITS. Additionally, we propose 
a slightly modified QACITS using two distinct image areas in¬ 
dependently. In Sect.[5l we report on simulation results of the 
QACITS performance in presence of higher order aberrations 
affecting the wavefront. In the flnal section, we draw the conclu¬ 
sions of our study. 

2. QACITS: Quadrant Analysis of Coronagraphic 
Images for Tip-tilt Sensing 

In this section, we introduce the QACITS post-coronagraphic 
technique to retrieve the tip-tilt affecting t he beam upstream a 
vortex phase mask. The demonstration by [Mas et al] (12012l) in 
the case of the FQPM is based on simulations and experimental 
data but an analytical model could also be derived (P. Baudoz, 



Fig. 2. Simulated images obtained for a tip-tilt of 0.2AID applied in the 
horizontal direction, for the case of a four-quadrant phase mask (left) 
and a vortex phase mask (right). Each quadrant Qt is a square of width 2 
AjD and defines an area where the flux is integrated in order to quantify 
the asymmetry in the image. 


private communication). In the present section, we derive the an¬ 
alytical model for the VC of charge L = 2, based on the typical 
coronagraph layout illustrated in Fig.tT] For that purpose, we use 
the Zernike formalism described in detail in Appendix lAl 

2.1. The quadrant analysis principle 

[Mas et al] (l2012h have shown that the amount of tip-tilt aberra¬ 
tion that affects the wavefront upstream the coronagraphic mask 
can be retrieved by analysing the residuals of the attenuated 
on-axis image acquired by the scientific detector. Indeed, this 
aberration induces an asymmetry in the pattern, as illustrated in 
Fig.|2l The asymmetry is quantified by two flux measurements, 
AIx and A/y, corresponding to the flux gradient along two or¬ 
thogonal directions in the image, which can be defined as 


,, (/2+/4)-(/l+/3) ,,, 

AIx = -and Aly 

lo 


(/l+/2)-(/3+/4) 

k 


( 1 ) 


with Ij^ = I the flux contained in each quadrant area Qk and 

4 . 

^ = Z n he the total amount of flux contained in the non coro- 

i 

nagraphic image /nc- In practice, these ar eas are squares o f width 
a few AjD {lAjD in Fig.O. Empirically. [Mas et"^ (12012h found 
that in the small aberration approximation, these quantities are 
directly linked to the amount of tip-tilt in the x and y directions, 
Tx and Ty respectively, following the model 


A4 = p pi + aT^p) and My = /3 {p + aTyT^, (2) 

P being a normalization coefficient. In the case of the FQPM, 
they found (Tfqpm = 4 and managed to find an approximate solu¬ 
tion of the system. In the following section, we show that for the 
vortex phase mask, the model can be derived analytically using 
the Zernike-based analysis. 


2.2. Zernike formalism: from entrance pupil to Lyot plane 

The phase of a tilted wavefront using Zernike polynomials ex¬ 
presses as 


0 = r,Z2(r) + r,Z3(r), (3) 

with r = (r, 0) the polar coordinates in the pupil plane, Tx and Ty 
the root mean square (rms) values for tip and tilt modes in radi¬ 
ans, Z 2 = 2r cos 6 and Z 3 = 2r sin 6 the tip-tilt mod es expmssed 
as the standard Zernike polynomials described by iNollI (Il976l) 
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and recalled in the Appendix lA.ll These polynomials are nor¬ 
malized to 1 radrms. In the small aberration approximation, the 
exponential function describing the wavefront can be expanded 
and approximated by 


F - 

^pup ■“ ^ 


0 ^ 


(4) 


The second order expansion is needed in order to make the non- 
symmetrical terms appear in the final PSF expression. The devel¬ 
opment of (p^ is a combination of the following terms, projected 
on the Zernike basis: 



Z2(r)2 = Zi(r) + Z4(r)/V3 + 2Z6(r)/V6, 

Ziirf = Zi(r) + Z4(r)/V3-2Z6(r)/V6, ( 5 ) 

Z2(r)Z3(r) = 2Z5(r)/V6, 

where Zi(r), Z 4 (r), Zs(y) and Z6(r) correspond respectively to 
the p iston, focus and the two astigmatism modes (see Appendix 
lA.ll for the explanation about the numbering of the polynomials). 
The complete expression of the field at the entrance pupil can 
thus be written as a linear combination of Zernike polynomials: 


-pup 


Ti + Ti 


Zi(r) -r iT^Z 2 {Y) + iTyZ^ir) 


^ X 

2 V 3 


-Z4(r) ■ 


2T,Ty 


V6 


ZsCr)■ 


Ti - Ti 


V6 


-Z 6 (r). 


( 6 ) 


As shown in Appendix IA.2I and IA.3I the decomposition of 
the wavefront onto the Zernike polynomial basis under the small 
aberration approximation turns out to be very convenient for de¬ 
scribing the effect of the vortex phase mask. Indeed, when prop¬ 
agating through the vortex focal plane mask to the Lyot plane, 
these Zernike modes translate into complex linear combinations 
of Zernike polynomials inside the geometrical pupil (the com¬ 
ponents outside the pupil are discarded here, as they are blocked 
by the Lyot stop). 

Using the conversion table that gives the field inside the re¬ 
imaged geometrical pupil after a VC of charge Ip = 2 (Table 
lA.lb . we can thus directly express the electric field after the Lyot 
stop as 


^Lyot — 


iTr 




-Z 2 (r) + 


■ iT^ 


-Z 5 (r) • 


Ti + Ti 


-Z 3 (r) - 


(T, + iTyf 

2V3 


2 V 6 


-Z6(r). 


■Z 4 (r) 


(7) 


2.3. Final image analysis 

The electric field in the detector plane is obtained by the Fourier 
transform of Eq.|7](the Fourier transform of the Zernike polyno¬ 
mials is reminded in Eq. lA.3l) . which leads to 


Z7 _ ^TT' I :'T' \2 2/3(27r(z) 

fidet - -(fx + U^) 

+ niT, + iTy)^-^^e'^ ( 8 ) 

1 2 /y 2 1 2 /^^^ 


Fig. 3. Simulation of the different components that form the final image 
on the detector when the input wavefront is tilted, with the notations 
used in Eq.|9] 


with (cr, if/) = a the polar coordinates in the image plane. The 
final image measured by the detector is the squared modulus of 
the previous expression and is thus given by 


/de, = [{tI + T^)A2{af + i [tI + Tff A^iaf 

+ ((r^ - Tp cos(2tA) + 2{TlTy + T^T,) sin(2^)) A^Ca) 

2 [iJl + T,T^) cosW + (T^y + TlTy) sinW) A 2 (a)A 3 (a)], 


(9) 


with Ai(a) = image consists of several terms, 

but only the last two contribute to the axial asymmetry in the 
image, as illustrated in Fig. [S] As a consequence, the relation 
between tip-tilt and the asymmetry simply writes 


A4=^(r3 + r,7'2), 

AIy=i3{T^y+TyT^). 

Here, JS is sl normalization constant, corresponding to 


( 10 ) 


P 


~h 



A2(a)A3(a)ada, 


( 11 ) 


with ao the maximal value allowed for a when integrating over 
the finite size quadrants (typically of few A/D). For the sake of 
simplicity, the integral has been expressed using polar coordi¬ 
nates, but rigorously it should be rewritten in order to consider 
the square shape of the quadrants. 

This result is fu l ly con sistent with the model found empir¬ 
ically by iMas et all (l2012h for the FQPM. The only difference 
lies in the value of the a parameter of the model given by Eq. O 
tJ^FQPM = 4 while (Tvortex = L In the vortex case, the system ad¬ 
mits a unique solution, which can be written as 


1 

1 A/2 


Ia/2 +A/2 

A/A5 

A/2 

P 1 

i A/2 + A/2 


( 12 ) 


The particular value of Ofyortex reflects the fact that a vortex phase 
mask is perfectly centro-symmetric, which is not the case for the 
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Fig. 4. Experimental results for the estimation of the tip-tilt aberration 
in one direction (the AGPM has been translated along the x axis only). 
Error bars are computed from the standard deviation of 50 values esti¬ 
mated from a sequence of 50 images. The images on the top left corner 
of the graph show the mean images of several sequences, acquired for 
different values of tip-tilt that are given on top of each image in A/D. 


FQPM. Indeed, from this system of equations, it can be shown 
that the simple law 

AIg=fiT^g (13) 

is sufficient to describe the relation between the tip-tilt and the 
asymmetry, as soon as the differential intensity AIo is measured 
along the axis of the applied tip-tilt Tq. In this case, the differ¬ 
ential intensity measured in the orthogonal direction, Ale+nji, is 
indeed zero. The direction of the tip-tilt can be inferred from 
the AIx and Aly measurements: tan^ = Aly/AIx, implying that 
AIq = (A/^ -r A/^)^/^. In other words, it means that the cross¬ 
terms TxTy and TyT^ in Eq. [T0|are not due to a cross-talk be¬ 
tween the two axes, like in the case of the FQPM, but are rather 
due to a change of reference axes. 

3. Experimental validation 

The model describing the relation between the tip-tilt amount 
and the asymmetry in the VC image has been validated thanks 
to experimental data. Test campaigns were indeed carried out for 
characterizing new-generation L-band AGPMs recently manu¬ 
factured at Uppsala Universitet (Vargas Catalan et al., in prep). 

These campaigns have been conducted on the YACADIRE 
bench at LESIA (Observatoire de Paris). This bench was 
used to characterize the coronagraphic masks for SPHERE 
(iBoccaletti et al.|[2Q0^ and thus mimics its optical layout (//40 
converging beam at the focal plane). We used a circular non- 
obstructed pupil and a circular Lyot stop (radius downsized 
by 80% with respect to the entrance pupil radius). For the 
testing of the AGPMs, a cold L-band spectral filter was in¬ 
stalled in the camera enclosure. The source is a Tungsten lamp, 
feeding a single-mocl e fibre. The bench layout is detailed in 
iDelacroix et al.l (l2Q13l) . who report on the first laboratory char¬ 
acterization of L-band AGPMs. 

The AGPM was first centred in x and y with respect to the 
beam by minimizing the flux integrated by the camera. The posi¬ 
tion along the optical axis was also optimized. Sets of 50 images 
were taken for different positions of the AGPM along the x axis. 
The AIx and Aly values are measured for every image. One has 


to note that translating the AGPM in the focal plane has not the 
same effect on the coronagraphic image as a tilted wavefront hit¬ 
ting the mask. Its shape will be affected in the same way, but it 
remains centred on the same position, while a tilted wavefront in¬ 
duces an additional translation of the image. This has been taken 
into account in the data processing (the quadrants were shifted 
by the number of pixels expected for the corresponding tip-tilt). 
The Tx and Ty are estimated for each image using Eq. [121 For 
one position, the final tip-tilt estimates result from the mean of 
the 50 estimates, and the error bar from their standard deviation. 

The results are shown in Fig. [4l The Tx estimates are in 
agreement with the true tip-tilt for a range of around +0.5 A/D 
from the center, where the estimations start to diverge from the 
expected value by more than their error bar. While Ty = Q was 
expected for the other axis, it seems that the position in that di¬ 
rection was not optimal and that the AGPM was probably off by 
about -0.07 A/D, corresponding to a shift of 10 [im in the focal 
plane. The data set has also been processed to estimate the trans¬ 
mission efficiency as a function of tip-tilt alon g the x direction. 
These results are detailed in the Appendix lB .2l and show that the 
highest extinction rate was obtained at 0.02 A/D from the posi¬ 
tion that was thought to be optimal during the experiment. This 
corresponds to a shift of 3.5 [im in the focal plane. 

In conclusion, our results show that the model derived to re¬ 
trieve the tip-tilt is valid for a circular non-obstructed pupil. The 
post-processing of the images has also shown that the manual 
optimization of the x and y position of the AGPM might not be 
optimal (for this particular experiment, the best manual align¬ 
ment of the AGPM was off by 0.02/I/D and 0.07/I/D in x and 
y, respectively), showing the limit of a manual positioning, as it 
is currently performed at the telescope. An automated method of 
tip-tilt retrieval based on the QACITS post-coronagraphic anal¬ 
ysis will thus significantly improve the vortex phase mask cen¬ 
tering. 


4. QACITS on a centrally obstructed pupil 

All the considerations from the previous sections are valid for 
a circular non-obstructed pupil. However, ground-based tele¬ 
scopes are usually centrally obstructed by the shadow of the 
secondary mirror. In the case of a central obstruction, the field 
distribution at the Lyot plane is affected by an additional contri¬ 
bution that falls inside the geometrical pupil, even for an on-axis 
source, thus preventing from a perfect on-axis starlight rejection. 
This significantly impacts the final image shape, and in partic¬ 
ular the asymmetry of the image. As illustrated in Fig. [5l the 
flux gradient changes sign for small tip-tilt in comparison with 
the image produced by an unobstructed pupil, implying that the 
model linking the differential intensity and the tip-tilt is differ¬ 
ent and more complex. In the following section, we analyse the 
theoretical model and adapt our QACITS tip-tilt estimator. 


4.1. Analytical derivation of the model 

Following the superposition principle, the entrance pupil can be 
written down as the sum of a positive contribution for the circular 
non-obstructed pupil and a negative contribution for the central 
obstruction of radius r < 1 (the pupil is defined with a radius of 
1 when using Zernike polynomials), that is 


r _ rO I robsc 

^pup “ ^pup ^pup 


0^00 _ 0^0obsc 


(14) 


Article number, page 4 of [T2| 

















E. Huby et aL: The QACITS technique 


a) Simulated images for a circular non obstructed pupil 



b) Simulated images for a circular obstructed pupil (24%) 


O 



Fig. 5. Each row of images corresponds to simulated coronagraphic 
PSFs obtained for different tip-tilt values, from left to right: 0.01, 0.05, 
0.10, 0.20 and 0.40 T/D. Each image intensity has been normalized by 
its maximum value. 


with 00 = 0 as defined in Eq. [S] The phase term 


^obsc — ^ 




(15) 


defines the phase of the component inside the central obstruc¬ 
tion, using scaled Zernike polynomials (defined for r/r < 1), 
such that the total field in the central obstructed area is cancelled 
out. As a consequence, the electric field at the Lyot plane (after 
the Lyot stop) will be composed of all the terms already men¬ 
tioned in Eq.[7]and of the additional following terms arising from 
the presence of the central obstruction: 


robsc 

^Lyot 


^2ie 


= + T^) - l) ^ {Ty + iT^) ^ 

J.2 , 




7^2 


, e4^‘^ 

-iT.Ty\—. 


3ie 


(16) 


Basically, these terms correspond to the decaying exponential 
term s that appear outside the ge ometrical pupil (see Appendix 
IA.3I and in particular Eq. lA.lSb of radius p, since the compo¬ 
nents inside the obstruction are blocked by the Lyot stop. As a 
consequence, these terms are defined for r > r and r < 1 (inner 
and outer diameter of the Lyot stop). The Lourier transform of 
a function of the general form restrained to this domain 

can be written as 


0.2 

.td 0.1 
c 

ro 0.0 
<- 0.1 
- 0.2 


- 3 - 2-10 1 2 3 

A/D 

Fig. 6. Horizontal profile (0 = 0) for the two asymmetric contributions 
in the final image plane, due to the whole pupil (i.e. A 2 {a)A 2 (a)) and to 
the central obstruction (i.e. AAi(a)A 2 (a)). 

1 - t(T^ -b Ty) ^ 1, thus assuming that the tip-tilt has a negligi¬ 
ble effect on the light diffracted by the central obstruction. As a 
consequence, the intensity on the detector consists of the expres¬ 
sion given in Eq.0augmented by the following terms (calculated 
from the modulus of the first term of Eq.[T8]and cross-terms be¬ 
tween the terms of Eq.[8]and first term of Eq.fT^: 

^ AAi(a)^ 

+ (T^, + T^)xAAi(a)A3(a) 

+ ((r^ - r2) cos(2^) + 2T,Ty sin(2^)) x AAi {a)A 3 {a) 

+ {27 cos(i//) + 2Ty sin((/>^)) X AAi(a)A 2 (a) j. 

(19) 

Only the last term of Eq.[T9]produces an asymmetric pattern with 
respect to the x and y axes. It has to be noted that the principal 
lobe of AAi(a)A 2 (a) has negative values, so that this term will 
be in competition with the asymmetric term arising from the cir¬ 
cular unobstructed pupil (last term of Eq.O. This is illustrated in 
Lig. [6l showing the horizontal profiles for each contribution. In 
addition, the contribution of the central obstruction is weighted 
by a coefficient directly proportional to the amount of tip-tilt, 
while the contribution of the circular pupil is lessen by the cube 
of the amount of tip-tilt. This explains why, for very small tip- 
tilt, the asymmetry in the images simulated with an annular pupil 
appears with a gradient of opposite sign compared with images 
simulated for an unobstructed pupil (Lig.O. Therefore, the rela¬ 
tion between the tip-tilt and the asymmetry in the image can be 
written as 



r 


^ikO 


= [Ak-i(a) - Ak-i(aT)\, 


(17) 


SO that the electric field on the detector due to the central ob¬ 
struction can be written as 


= ^(l+ 7’2))e2''^ r^AAiCa) 

+ n{iTy-T,)t^‘^ T^AA2{a) (18) 

+ T^AA^ia), 

with AAkia) = Ak{a)-Ak{aT)lT^~^. Numerical estimations show 
that the AA\t^ component is the dominant term, the two other 
ones being significantly smaller in absolute values due to the 
factor (because r < 1). Since it would be uselessly painful to 
derive the complete expression of the intensity recorded by the 
detector, we choose to neglect the two weaker terms in the fol¬ 
lowing computation. In addition, we can approximate the factor 


My=l3{Tl + TyTl) + yTy, 

with p and y two real parameters of opposite signs. As it will 
be illustrated in the following section, the main issue with this 
model is that it will necessarily limit the range where the stan¬ 
dard QACITS method can be applied, because the competition 
between the two terms leads to a possible ambiguity to retrieve 
the tip-tilt from a single intensity measurement. It also reduces 
the sensitivity, as the two contributions partially cancel each 
other. That is why a dual area QACITS method is proposed in 
the next section. 

4.2. QACITS in dual areas 

As illustrated in Lig. [ 6 l the asymmetric contribution due to the 
central obstruction undergoes a sign inversion at 1.6 AjD, which 
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b) Flux repartition between the 
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Fig. 7. a) Simulated image for a eentrally obstrueted pupil (24% of the 
diameter pupil) with 0.2 A/D tip-tilt, showing the spatial flux distribu¬ 
tion for the standard, inner and outer areas. The total image width is 
6 A/D while the circle has a radius of 2 A/D. h) The flux repartition be¬ 
tween the inner (blue eontinuous line) and the outer (red dashed line) 
areas as a function of tip-tilt (estimated from simulated images). 


corresponds to the zero of the A 2 (a) function. Two areas can 
thus be defined in the image: the inner area (r < 1.6 A/D) and 
the outer area (r > 1.6 A/D). Note that when the outer diame¬ 
ter of the Lyot stop is downsized, this 1.6 A/D boundary has to 
be scaled proportionally (for instance, a Lyot stop downsized by 
80% has for effect to push the boundary to 2 A/D). These ar¬ 
eas are shown on a simulated image in Fig. [7^, highlighting the 
fact that the intensity gradient has opposite sign depending on 
the considered region. Fig.fTJ) shows the flux repartition between 
these areas for an unobstructed circular pupil and an obstructed 
pupil. While in the ideal unobstructed case, the flux is mostly 
concentrated in the central lobe (~ 80% of the total flux), a sig¬ 
nificant portion of the flux spreads outwards in presence of a 
central obstruction, making this dual measurement legit. 

The differential intensities corresponding to the standard 
QACITS, and to the QACITS split down into inner and outer 
areas are shown in Fig. [8] for the case of an annular pupil (24% 
obstruction in diameter). While the differential intensities com¬ 
puted in the standard way show a degeneracy and a limited am¬ 
plitude, the intensities restricted to the inner and outer areas 
reach higher absolute values. An interesting feature appearing in 
these plots is the fact that for small tip-tilt (< 0.2A/D), the model 
can be approximated by the linear part of the model, which dom¬ 
inates over the cubic term. For the sake of simplicity, we will 
use this approximation thereafter, especially since the system of 
equations (Eq.[20|) does not lead to simple analytical solutions. 

The p and y parameters defining the model given in Eq. [20| 
have been estimated thanks to simulations for different pupil 
configurations, and in particular different Lyot stop parameter 
values: the inner and outer diameter, Lin and Lout, defined as a 
fraction of the entrance pupil diameter D. The p and y parame¬ 
ters correspond to the cubic and linear components respectively 
and are computed by fitting the simulated points in the least- 
squares sense. The values are reported in Table [U These results 
show that the y parameter weighing the linear part of the model 
increases with the reduction of the outer diameter of the Lyot 
stop mask. This is expected since the flux due to the diffraction 
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I 0.00 


- 0.01 

- 0.02 

-0.03 


0.0 0.1 0.2 0.3 0.4 

Tip-tilt (A/D) 

Fig. 8. Estimated differential intensities resulting from simulated im¬ 
ages in the case of a pupil centrally obstructed (24% of the full diam¬ 
eter). The solid line curves show the best fit model in the least squares 
sense (the model consists of a linear and a cubic component), while the 
dashed lines show the linear contribution only. The three cases differ in 
the area used to compute the differential intensity: standard whole area 
(in black), inner (< 2A/D, in blue) or outer area (> 2A/D, in red). The 
outer diameter of the Lyot stop is downsized by a factor of 80%, while 
the central obstruction diameter is set to 35% (1.45 oversizing factor). 

Table 1. p and y parameters estimated from simulations and based on 
the whole central image (standard method) or only on the inner or outer 
areas of the image. Different Lyot stop configurations have been simu¬ 
lated, with the first line corresponding to an unobstructed entrance pupil 
(hence the non specified inner diameter of the Lyot stop), while all the 
other cases result from an annular entrance pupil with 24% central ob¬ 
struction. 
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-^n 
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In. area 

Out. area 

(%) 
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p 

r 

Pin 

Tin 

Poxxi 

Tout 

100 

/ 

1.08 


0.94 


0.07 


100 

24 

0.93 

-0.04 

0.88 

-0.10 

0.04 

0.06 

100 

35 

0.90 

-0.04 

0.79 

-0.10 

0.11 

0.05 

80 

24 

0.75 

-0.05 

0.70 

-0.13 

0.05 

0.08 

80 

35 

0.68 

-0.05 

0.56 

-0.11 

0.12 

0.06 


by the central obstruction mainly distributes to the area close to 
the central obstruction (see the decaying exponential functions 
of Eq.O, while the tip-tilt energy coming from the whole pupil 
is spread over the whole pupil. As a consequence, cropping part 
of the outer rim of the pupil implies that the central obstruction 
contribution, which is the source of the linear dependency, be¬ 
comes relatively stronger. 

The parameters have also been estimated for measurements 
restricted to the inner and outer areas. In both cases, the y param¬ 
eter reaches higher values, and thus provides a better dynamic, 
in comparison with the values obtained by integrating the flux in 
the whole image (standard method). The final estimator is there¬ 
fore taken as the average of the inner and outer estimators based 
on the linear approximation of the model, and can thus be written 
as 


2 \ Tin Tout / 


( 21 ) 


with and A/ defined as vectors with the v and y components 
of the tip-tilt estimate and differential intensity measurements re¬ 
spectively. The exponents "in" and "out" refer to the area of the 
image used to integrate the flux, namely inner or outer part. This 
average estimator has been applied to simulated images affected 
by a tip-tilt ranging from 0 to 0.4 A/D. The tip-tilt residuals are 
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Tip-tilt (A/D) 

Fig. 9. Tip-tilt residuals obtained for the same simulation parameters as 
Fig.ia.e. 24% central obstruction, Lyot stop of 35% and 80% inner 
and outer diameter). 

reported in Fig. [9l The mismatch between the model and the 
linear approximation induces a bias in the inner and outer es¬ 
timators. These biases happen to be of opposite sign, and thus 
compensate each other at least partially when taking the aver¬ 
age. In practice, this offset is not critical since the QACITS al¬ 
gorithm is supposed to be used in closed loop control. The re¬ 
sults show that for large tip-tilt amounts, the combined estima¬ 
tor under-estimates the amplitude, which means that the conver¬ 
gence might be slower at first. 

To conclude, we have derived the theoretical model and mod¬ 
ified the QACITS estimator to make it applicable to the case of a 
centrally obstructed aperture. Because the contribution of the ob¬ 
struction counterbalances the contribution of the circular pupil, 
the validity range is reduced to small tip-tilt amounts (for tip-tilt 
< 0.2/l/Z), the bias is smaller than 3%). However, the presence of 
the central obstruction is responsible for a higher starlight leak¬ 
age (at least 5% for a central obstruction of 24% in diameter), 
providing a better sensitivity, but also a better dynamic due to 
the linearity of the model, as opposed to the cubic model in the 
non obstructed case. 


5. Performance in presence of higher order 
aberrations 

In practice, real wavefronts are affected not only by tip-tilt but 
also by higher order aberrations. Static aberrations due to im¬ 
perfect optics surfaces can be handled by subtracting a reference 
image. However, quasi-static speckle patterns may corrupt the 
AIx and Aly measurements. Such aberrations may be caused by 
temperature and mechanical drifts, that slowly evolve with time, 
and are not sensed by the adaptive optics system. In order to 
quantify the effect of higher order aberrations, simulations were 
conducted with phase screens generated from a power spectral 
density defined as the inverse power law of exponent 2. This 
kind o f model is typical for fractal finish surface quality (IChurchl 
1 19881). like the high quality optics of the SPHERE instrument 
( Dohlen et al .1120lib . The simulated coronagraph is based on a 
circular entrance pupil obstructed by 24% of its diameter and a 
Lyot stop with an oversized central obstruction of 35% and outer 
diameter of 80% of the initial entrance pupil (this corresponds to 
the typical obstruction and Lyot configuration of the NIRC2 in¬ 
strument at the Keck telescope). 

Every phase screen is drawn randomly. The tip-tilt compo¬ 
nent is estimated by a projection onto the base of Zernike poly¬ 
nomials and subtracted. A hundred tip-tilt values ranging from 0 


to 0.4 T/D are uniformly drawn and applied to the wavefront in 
the horizontal direction (orientation angle ^ = 0). The tip-tilt am¬ 
plitude and orientation angle are then estimated from the image 
with the dual QACITS method. The final estimate is computed 
from the average of both estimators using the inner (<2 A/D) and 
the outer part of the image (2 A/D<a<3 A/D). As shown in the 
previous section, for small tip-tilt values (<0.2 A/D), the relation 
between the asymmetry in the image and the tip-tilt amount can 
be approximated by a linear function, whose proportionality fac¬ 
tors, 7 , are given in Table [T] i.e. = -0.11 and yout = 0-06 in 
the configuration used in our simulations. 

The results are shown in Eig.[T0l Different aberration levels 
have been simulated, from Srms = T/IO"^ to Srms = T/IO. As 
expected, for large tip-tilt amounts (>0.2 A/D), the model is not 
valid any more and a bias appears, in the amplitude estimation in 
particular. The root mean square (rms) values of the tip-tilt am¬ 
plitude residuals reported in Eig. [TT]have therefore been com¬ 
puted on the reduced range of tip-tilt < 0.2A/D. These results 
show that in the small tip-tilt regime and for very low aberra¬ 
tion levels (Srms < A/300), the bias due to the linear approxi¬ 
mation dominates the speckle noise, and limits the accuracy of 
the estimation to 2.2 x 10“^ A/D. Eor higher aberration levels, 
the accuracy is dominated by the effect of the aberration, and the 
tip-tilt residual rms increases linearly with the wavefront error 
rms, expressed as a fraction of wavelength. This is observed for 
the amplitude as well as for the orientation angle (the slope of 
the best fit models drawn in Eig. [TT]in log-log scale is 1.0 for 
both cases). 

These results illustrate the stability level we can expect from 
a control loop based on the QACITS technique when higher or¬ 
der aberrations affect the PSE shape. The tip-tilt affecting the 
beam can be estimated with a precision better than 10“^ A/D and 
5 X 10“^ A/D in presence of wavefront errors up to ^rms = A/10 
and ^rms = A/14 respectively (corresponding to ~ 50 nm rms and 
-270 nm rms at 3.75[im). However, quasi-static speckles tend to 
evolve slowly with time (i.e. on minute time scales). Therefore, 
in practice, two consecutive images are not completely decorre- 
lated, unlike our set of simulated phase screens, and in this case 
part of the high order aberration impact can be avoided by sub¬ 
tracting a reference image, obtained for the best centering of the 
coronagraphic mask. 

6. Conclusions and prospects 

We have described the QACITS technique for the vortex coro¬ 
nagraph, a method ori ginally introduce d in the case of the four- 
quadrant phase mask (iMas et aH l2012h and a circular non ob¬ 
structed aperture. We have derived the analytical model for the 
VC and found a cubic power law, validated by simulations and 
experimental results. However, the presence of a central circu¬ 
lar obstruction adds a linear component that induces an intensity 
gradient in the opposite direction. In order to tackle this more 
complex model, we have introduced the QACITS method in dual 
zones (distinguishing the inner lobe from the external region), 
which allows the disentanglement of the cubic and linear com¬ 
ponents. 

Simulations of a typical telescope configuration carried out 
in the presence of higher order aberrations show that the QAC¬ 
ITS method provides an estimation of the tip-tilt with a precision 
of 10~^ A/D for wavefront errors amounting to A/70 rms. Eor 
very low level of aberrations (< A/300), systematic errors arising 
from the linear approximation of the model limit the accuracy of 
the estimation to 2.2 x 10“^ A/D. The practical implementation 
may also be limited by other factors, such as the brightness of 
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Fig. 10. Simulation results of tip-tilt estimation using QACITS with a centrally obstructed pupil in presence of higher order aberrations. Amplitude 
residuals are shown on the left, while orientation angle residuals are shown on the right, for different levels of aberrations. 
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Fig. 11. Root mean square values for the residual tip-tilt amplitude (left) and orientation angle (right), as a function of the wavefront error. Note 
that for the amplitude, the rms is computed over the reduced [0,0.2]/l/D range, where the linear model approximation is valid. The green dashed 
lines correspond to a best fit power law model. 


the star, or the possible asymmetry of the observed object. This 
aspect will be discussed in more detail in another paper. 

It can also be emphasized that the Zemike-based analysis 
reported in the Appendices highlights a remarkable feature of the 
vortex coronagraph: at first order, small aberrations expressed 
as Zernike polynomials simply translate into a complex linear 
combination of other Zernike polynomials in the Lyot plane. We 
are currently investigating other wavefront sensing techniques 
exploiting this characteristic. 

To conclude, the QACITS technique offers an easy way to 
control the centering of the vortex phase mask directly from the 
scientific image, thus avoiding non-common path errors that an 
additional wavefront sensor fails to measure. Its simplicity of 
implementation makes QACITS a valuable and directly available 
tool for all the instruments equipped with a vortex phase mask. 
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Appendix A: A Zernike-based analysis 

We propose a Fourier-based analysis of beam propagation using 
Zernike polynomial decomposition of the wavefront. Our com¬ 
putations are based on the standard layout of a coronagraph, il¬ 
lustrated in Fig. [U 

Appendix A.1: The Zernike polynomials 

The Zernike polynomials were described bv iNollI (Il976h and are 
defined for r < 1 (r = (r, 6) are the polar coordinates) as 


At the focal plane, the vortex phase mask ind uces a phase 
shift d epending on the azimuthal angle if/. Indeed, iMawet et~^ 
(I 2 OO 5 I) have shown that for a perfect vortex phase of topologi¬ 
cal charge Ip, the right- and left-handed circular polarization unit 
vectors are translated into left- and right-handed circular polar¬ 
ization vectors respectively, and are affected by a phase ramp 
^iipif/ Q-iipif/ respectively. The coronagraphic effect will occur 
for any value of Ip that is even. 

To ease the comparison between the field in entrance pupil 
(E’pup) and Lyot plane (E’Lyot), an inverse Fourier transform, noted 
is finally applied, leading to 


Zj{r) = y/^r+lR’^ir) y/2C„,{0) for m 0, 
Zj(r) = ^In + lf?°(r) for m = 0, 


with 


(n-m)/2 

Kir)= Z 

s=Q 


{-iy{n-s)\ 


5 '![(^ -I- m)/2 - 5 ’]![(^ - m)/2 - 


-r” 


-2s 


(A.2) 


Here, n and m are non negative integers satisfying m < n, with 
n - m even (in other words, n and m have the same parity). The 
azimuthal functions, CmiO) = cos mO and CmiO) = sin mO, are de¬ 
fined for even and odd j respectively, thus corresponding to real 
symmetric and antisymmetric modes respectively. The index j 
is a usual numbering system for the different modes, that will be 
used in the following developments. There is no equation linking 
the index j and the (n, m) integer pairs. The correspondences for 
the first 14 polynomials and their usual aberration designation 
can be found in Fig. IA.il 

The Fourier transform of the Zernike polynomials, noted Zj, 
can be written as 


Zj(a) = V2(^ + l)7rC^(0^)/ fQ^m 0, 

Ina 

Zj(a) = yn-\- l 7 r(-l )2 --- for m = 0, 

2na 

(A.3) 

with (a, ij/) = a the polar coordinates in the conjugate plane, and 
Jn{(T) the Bessel function of the first kind. 

Appendix A.2: The small aberration assumption 

Under the small aberration hypothesis, the wavefront at the en¬ 
trance pupil (amplitude of 1, phase 0) can be directly approxi¬ 
mated at first order as a combination of Zernike polynomials: 


^pup = exp(/0) ^ I+ i(p = Zi(r) + / Y^ajZjir), (A.4) 

j=2 

where r = (r, 6) are the polar coordinates in the pupil plane and 
a j a set of real coefficients descri bing the aberrations. 

The Fourier transform of Eq IA.4I leads to the field distribu¬ 
tion in the focal plane and thus consists of a linear combination 
of Zernike polynomial Fourier transforms, Zj. We can thus write 


^Lyot = +iZajr-^[Zj(aW‘>’'^] 

'f (A.6) 

= ^i(r) +/ Z 

j=2 

with denoting the field distribution in the 

Lyot plane when the input pupil amplitude is defined by the 
Zernike polynomial Zj. 

The first term fi(r) results from the perfect plane compo¬ 
nent (piston mode) that is completely diffracted outside the ge- 
ometric pupil in th e Lyot plane as long as the charge /p is even 
(iMawet et al.ll2005l) . In the following section, we derive the gen¬ 
eral expression of ^j(r) and show that they can be expressed as 
Zernike polynomials inside the geometrical pupil. 

Appendix A.3: The conversion tables Zj 

The expression of can be expanded using Eq. IA.31 thus be¬ 
coming 


6-(r) 


^ r- 2Jn,,(2na) 

" Jo 27 ro' Jo 


(A.7) 


with 


A’^= V2(«+ l)7r(-l)'^r'” form?t0, 
A® = V« + \n{-\)2 form = 0, 

and 


(A.8) 


Cm 


cos(m(A) = 

< sin(OT^) = f (e™'^ - 

1 


for m 9 ^ 0 and even j, 
for m 9 ^ 0 and odd j, 
for m = 0. 

(A.9) 


^j(r) will thus be written as one term (m = 0) or as the sum of 
two terms (m ^ 0). In any case, all these terms have the same 
form and can be simplified using the integral form of the Bessel 
function that expresses as 


1 

^ Jo 

which can also be written in the more convenient manner 


^foc ~ '^i(tr) + i ^ CLjZj{(x) 
J=2 


(A.5) 


p27T 

Jo ‘ 




2nfe’^^Jk{z), 


(A.11) 
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a) Charge Ip = 2 


b) Charge Ip = 4 



Fig. A.l. Qualitative representation of the field distribution at the Lyot plane for a VC of charge Ip = 2 (left) and Ip = 4 (right), as a function of the 
(n, m) integer pair of the input Zernike polynomial Zj. In the column m = 0, only one term defines the field distribution, and is depicted in green. 
For all the other cases, two terms contribute to the electric field, represented in cyan and pink. In the case where both terms are superimposed 
either outside or inside the pupil, these colors are mixed into a purple hue (see for instance the coma for the charge Ip - 2 vortex, or the tip-tilt for 
the charge /p = 4 vortex). 


where we identify z = 27ira, and k = Ip when m = 0 or k = lp±m 
when m ^ 0. 

Replacing the integral over the variable if/, the general ex¬ 
pression of ^jir) becomes 




2jh 


for m 9 ^ 0 and even 7 , 

for m 9 ^ 0 and odd 7 , (A. 12) 

for m = 0 , 


with 


= 



Jn+\(2na)Jk(27rar)da. 


(A.13) 


According to the Eq. 9 of iNollI (Il976l) . the integral can be 
linked to the R^(r) function (recalled in Eq. lA.2l) . that we rewrite 
here in a more general manner (not restricted to the conditions 
k <n and k and n with the same parity): 


for 1^1 < ^ j 


for 0 < r < 1 

10 

for r > 1 

1 

[0 

for 0 < r < 1 

for 1^1 > ^ j 


for r > 1 


(A. 14) 


One can note that this integral is discontinuous. The general form 
of ^j(r) can thus be considered as two components that are non¬ 
zero exclusively inside or outside the geometrical pupil, depend¬ 
ing on the comparison of |/^ + m| with n. We can thus write 


(A.15) 


where and f represent the contributions of the field inside 
and outside the pupil respe ctively. This field distribution is visu¬ 
ally illustrated in Eig. IA.il for a VC of charge Ip = 2 and Ip = 4. 
Note that or can be zero. Indeed, if m = 0, there is 
only one term (Eq. lA.12l) . which is defined either for r > 1 or 
0 < r < 1. In particular for the piston term (plane wavefront), 
there is no component inside the geometrical pupil, confirming 
the theoretical perfect extinction of the VC. Another interesting 
example is the case of the defocus and tip-tilt modes: in the case 
of the VC of charge Ip = 4, both of them fall outside the geo¬ 
metrical pupil, while there is a non zero contribution inside the 
pupil for the VC of charge Ip = 2. This result confirms that, for 
circular unobstructed pupils, charge 4 vortices are less sensitive 
to tip-tilt and defocus aberrations than charge 2 vortices, since at 
first order, these modes are completely rejected outside the pupil. 

The final expression of f/r) thus depends on the value of m 
and the parity of 7 . As a summary, we can write: 

- when m = 0 , it implies that n is even and hence: 




for |/p| < «, 0 < r < 1 , 

- V«+T 7 e'''’®/?j;+ij(l/r) for \lp\ > n, r > 1 . 

(A.16) 


“ when m 9 ^ 0 , we distinguish the cases even and odd 7 : 


for even 7 , 
+ for odd 7 . 

with 


(A.17) 




for 1^1 < ^, 0 < r < 1 , 


I aikO 


(A. 18) 


0 = + 
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Table A.l. Conversion table for the first eight Zernike polynomials for 
a charge Ip - 2 vortex phase mask. Only the contribution inside the 
geometrical pupil is considered (electric field after the Lyot stop). As a 
reminder, corresponds to the contribution of the input Zernike poly¬ 
nomial Zj inside the geometrical pupil, such that the table should be 
read line by line (for instance, if the entrance pupil contains the tip-tilt 


mode Z 2 , this will translate in the Lyot plane as ^2 

n _ 

Z 2/2 + iZilX). 


Z 2 

Z 3 

-^4 ^5 '^6 
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2 

1 
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2 

Table A. 2. Same as table lA.ll but for 
mask. 

a charge Ip 

= 4 vortex phase 
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Z 4 Z 5 

Ze 

Z 7 Zs 
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It is interesting to note that, if they exist, the terms that are 
defined inside the geometrical pupil (r < 1) can be expressed as 
a complex combination of Zernike polynomials. We can indeed 
write 


” ifk = lp±m = 0, 

(A. 19) 

with the + sign corresponding to the sign of k. As a consequence, 
a conversion table can be established, that gives the coefficients 
of the Zernike polynomials defining the field after the Lyot stop 
(i.e. since the is blocked by the aperture stop) for a given 
input Zernike polynomial, Z^, passing through the VC. These 
tables are given for the charge /p = 2 and Ip = 4 vortices (Table 
lAHandlA^. 


off-axis source close to the center. Because of the central sym¬ 
metry, only one axis is needed to describe the transmission as a 
function of the distance from the axis, noted T in rad rms. At the 
first order, the total transmission is estimated from Eq.[7]by 


mMT) = 


fup 


.yotl 


Xup 


T2 

T' 


(B.l) 


For convenience, the tip-tilt rms in radian, T, can be converted 
into an amplitude S in unit of AjD by means of the relation 
7^[rad] = s YxjD] X tt/Z, leading to 


m ,=2 (S) 



2 

WDY 


(B.2) 


This result is slightly different from the formula given by 
iJenkinsI (l2008l) . who derived it empirically. His result is also 
based on a square law but the multiplicative factor is different 
{n^/6 instead of tt^/S). Simulations have been performed in order 
to compare the two models. The parameters of the simulations 
are: a grid size of 1024 points in width, entrance pupil covering 
102 pixels and a Lyot stop of the same size as the entrance pupil. 
Particular care has to be given to numerical errors: the main part 
of them can be avoided by computing the entrance pupil profile 
that leads to perfect attenuation of an on-axis source. This is per¬ 
formed by simulating the propagation of a perfect circular wave- 
front up to the Lyot plane, cancelling out the residuals inside the 
geometrical pupil (relying on the argument that this is true ana¬ 
lytically), and final ly propagating th e result backwards, down to 
the entrance pupil (iKrist et al.ll2012h . The complex profile of the 
entrance pupil obtained in this way is used as the perfect wave- 
front. The results of the tip-tilt simulations are shown in Fig. IB. II 
and confirm that for small tip-tilt values, the transmission effi¬ 
ciency follows the function given in Eq. IB. 21 


Appendix B.2: Experimental results 

The experimental data described in Sect. [3] have been processed 
in order to estimate the transmission efficiency as a function of 
tip-tilt. The flux has been integrated for each position of tip-tilt in 
a square of width 10 A/D centred on the PSF, and divided by the 
value obtained for the AGPM translated by 7 A/D, a distance at 
which the beam is barely affected by the vortex phase mask. The 
transmission curve is shown in Fig. IB. 21 A polynomial function 
(composed only of even orders up to the 6^^, because of the ob¬ 
vious and expected symmetry) has been fitted to the data points. 
The best fit model leads to a position of the minimal transmission 
around -0.02 A/D, meaning that the position that was thought to 
be the optimal was actually off by 3.5 [im in the focal plane. The 
inner-working angle, defined as the distance where the off-axis 
transmission reaches 50%, is estimated to be 0.9 A/D (with D the 
diameter of the entrance pupil). The results were also compared 
to the theoretical model as derived in Eq. IB. 21 but the sampling 
was obviously not sufficient at very small tip-tilt to perform a 
useful comparison. 


Appendix B: Off-axis transmission 

Appendix B.1: Analytical function 

The Zernike analysis that has been carried for a tilted wavefront 
(Eq.U]) allows the estimation of the transmission efficiency for an 


Article number, page 11 of [12] 



















A&A proofs: manuscript no. Huby_2015_QACITS 



Tiptilt (^i/D) 


Fig. B.l. Transmission efficiency as a function of the angular distance 
from the center. The crosses result from simul ations while th e dashed 
line shows the theoretical function as stated bv iJenkinsI (l2008h . and the 
solid line shows the model derived in this work. The two are in agree¬ 
ment for the square dependency, but differ in the multiplicative factor, 
equalling ;r^/6 and tt^/S respectively. 



Tip-tilt (A/ D) 


Fig. B.2. Measured transmission efficiency for an off-axis source 
through a VC. The red circles correspond to experimental results, while 
the dark dashed line is the best fit model (polynomial function of even 
orders only, up to the 6th order). The li ght gr ay dashed line corresponds 
to the theoretical model as stated by Eq. IB. 21 assuming very small tip-tilt 
(it is thus drawn only for absolute tip-tilt < 03 A/D). The inner work¬ 
ing angle is graphically represented by the dotted lines, that highlight 
the transmission limit of 50%, reached for tip-tilt of 0.9 A/D. 
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